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Abstract 

Based on an effective chiral theory of pseudoscalar, vector, and axial-vector mesons, 
the coefficients of the chiral perturbation theory are predicted. There is no new param¬ 
eter in these predictions. The current quark masses, rri^. fx, and f v are determined 
too. 
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The chiral perturbation theory(ChPT) proposed a decade ago[l] is successful in parametriz¬ 
ing low-energy QCD. The chiral symmetry revealed from QCD , quark mass expansion, and 
momentum expansion are used to construct the Lagrangian of ChPT 

C = SgrrDpUD^ ± &Tr X (U ± U ] ) ± L 1 [Tr(T» M C/^C/ t )] 2 + L 2 (TrD ll U D V U^) 2 
±L 3 Tr( J D% J D / %) 2 ± L 4 Tr{D^UD^)Tr X (U ± U ] ) ± L 5 Tr £>%£>% { X U ] ± U X ) 
+L 6 [Tr X (U ± %] 2 + L 7 [Tr X (U - %] 2 ± L s Tr( X U X U ± X U\U^) 
-iL 9 Tr(F^D»UD u U' ± F*D' i U*D v U) ± L w Tr(F^UF^ R U ] ). (1) 

The parameters in the chiral perturbation theory(l) are pion decay constant f n and the 10 
coefficients. These parameters are determined by fitting experimental data[l,2,3,4,5,6,7]. 


Table 1: Values of the coefficients[1(b)] 
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The chiral perturbation theory is rigorous and phenomenologically successful in describing 
the physics of the pseudoscalar mesons at low energies. Models attempt to deal with the 
two main frustrations that the ChPT is limited to pseudoscalar mesons at low energy and 
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contains many coupling constants which must be measured. The price paid for including 


more mesons and determining the coefficients is to make additional assumption. Many 
models try to predict the chiral couplings(see Table 2 [8]). As pointed in Ref. [8], the most 
reliable one appears to be the influence of the low-lying resonances[4], There are many 
attempts to construct an effective theory containing vector and axial-vector mesons too[9]. 
Obviously, a successful extension should recover the results obtained by the ChPT and 
predict the coefficients of the ChPT. 


Table 2: Coefficients obtained by models 
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Nucleon loop 
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I have proposed an effective chiral theory of pseudoscalar, vector, and axial-vector mesons[10,ll]. 
It provides a unified study of meson physics at low energies. It is used to predict all the 
coefficients of the ChPT in this paper. Before doing that it is meaningful to present more 
detailed arguments about the ansatz made in this theory. The ansatz is that the meson fields 
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are simulated by quark operators. For example, 


P\ L = - (2) 

A 9pm 2 p 

The ansatz( 2 ) can be tested. Applying PCAC, current algebra, and this expression( 2 ) to 
the decay p —> 7 T 7 t, under the soft pion approximation it is derived 

2 fpmrdp = 1 (3) 


which is just the result of the VMD[12], As a matter of fact, using fermion operator to 
simulate a meson field has a long history. More than six decades ago, Jordan et al [13] 


observed 

- 4=3 P<t> = ^ 757 ^ 

V 71 


(4) 


in 1 + 1 field theory. A similar relation as Eq.(4) is proved in the bosonization of 1+1 
field theory[14], On the other hand, using quark operators to simulate meson fields have 
been already exploited in the Nambu-Jona-Lasinio(NJL) model[15], a model of four quark 
interactions. 

In Refs. [10,11] the simulations of the meson fields by quark operators are realized by 
a Lagrangian which is constructed by using chiral symmetry and the minimum coupling 
principle 


£ = ■ d + 7 • v + 7 • ay 5 — m«(a;))^(a;) — 
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( 5 ) 



where M is the quark mass matrix 


^ m,, 0 0 


0 rn d 0 


y 0 0 m s j 


v » ~ + KK* a + (§ + 7^3)07, + (| — -^A 8 )0^, — Tja^ + A a A^ + (§ + 7^3)/^ + (A 

73 A s and u = exp{i^(jiT\i + \ a K a + A 8 ? 7 8 + 770 )}. u can be written as 


U = ^(1 + 75 )u + i(l + 75 )^, 


( 6 ) 


where U = exp{i{ri^i + A a K a + X 8 p + 70 )}• Since pions, kaons and rj are Goldstone bosons 
we treat them differently from vector and axial-vector mesons. Through the scheme of 
the nonlinear er model the pseudoscalar mesons are introduced to the Lagrangian(5). This 
treatment is backed by the ChPT. In the Lagrangian of the ChPT(l) the first two terms 
are from the nonlinear a model. The matrix representations of the pseudoscalar mesons, U 
and U\ in Eq.(5) are the same as the ones in the ChPT(l). Spontaneous chiral symmetry 
breaking is revealed from this mechanism [10] when taking u 1. On the other hand, mesons 
are bound states solutions of QCD they are not independent degrees of freedom. Therefore, 
there are no kinetic terms for meson fields. The kinetic terms of meson fields are generated 
from quark loops[10,ll]. Using the least action principle, the relationship between meson 
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fields and quark operators are found from the Lagrangian(5). Taking the case of two flavors 
as an example, from the least action principle 


dC dC dC dC dC dC 

— = 0 , — = 0 , -- = 0 , — = 0 , - = 0 , — = 0 , 

ollj dp dp^ da \i o/m 

we obtain 


(7) 


— = 75^ + ixi\)Tii>) / 

a 

x = (i' 075'0 - — ipTiip)/ (ipijj + i—^ 75 ^), 

O' O' 

Pm = < = --^wy 5 ^, 

m 0 m 0 

1 - 1 

=-/> =- 2 ^ 7 ^ 75 - 0 , ( 8 ) 

7715 mo 

where a + iysT ■ fl = ue~ trn% , o' = a/ 1 — If 2 , and x = tam 7 . The pseudoscalar holds have 
very complicated quark structures. Substituting the expressions( 8 ) into the Lagrangian(5) 
of two flavors, a Lagrangian of quarks is obtained, ft is no longer a theory of four quark 
interactions. The vector and axial-vector parts are the same as the ones in the NJL model. 
However, the pseudoscalar part is different from the NJL model. 

As shown in Refs.[10,11] there are explicit chiral symmetry(PCAC[ 19] ), dynamical chiral 
symmetry breaking, and large Nc expansion. They originate in QCD. However, theoret¬ 
ically, we do not know how to prove the ansatz of simulating the meson fields by quark 
operators for 1 + 3 QCD. We use the Lagrangian(5) to study meson physics at low en¬ 
ergies to see whether the Lagrangian is phenomenologically successful. In Refs. [10,11] the 
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Lagrangian(5) has been used to calculate the masses and electromagnetic, weak, and strong 
properties of the pseudoscalar, vector, and axial-vector mesons. Theoretical results agree 
well with the data. Taking the masses of the p and the a\ meson as an example of the 
results obtained by this theory. The mass of the p meson is originated in dynamical chiral 
symmetry breaking[16]. The mass difference of the p and the a\ mesons is from spontaneous 
chiral symmetry breaking[10]. 

Large Nq expansion is part of this effective chiral theory[10]. All the calculations are 
done at the tree level and meson loops are at higher order in large Nq expansion[10]. All the 
masses of the mesons are less than the cut-off[10]. Therefore, this theory is self-consistent. 

By using the ansatz, this effective chiral theory extends the ChPT to include vector and 
axial-vector mesons. However, this theory faces the second challenge: predicting all the 
coefficients of the ChPT. This is the task of this paper. The study done in this paper is at 
the tree level. 

As a matter of fact, in Ref.[10] two of the coefficients have been determined from the 
study of 7T7T scattering 


CX 2 — —— C , 


C = 


f 

J 7T 


2 gml' 


(9) 

( 10 ) 


where «]. 2 are defined in Ref. [4] and g is the universal coupling constant of this theory[10,ll]. 
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Comparing with Eq.(l), we obtain 


OL 1 — ALi + 2E3, «2 — AL/ 2 - (11) 

Eq. (9) is the same as the one obtained in Ref. [4] in which the p resonance participates in nn 
scattering. From Eq.(ll) it is obtained 

2(Li + L 2 ) + L 3 = 0, (12) 

U = p 2 - (13) 

The decay A ")4 has been used to determine the coefficients of ChPT[17,18]. Both vector 
and axial-vector currents contribute to K\\. In this paper we only try to predict A 3 from 
Ki 4 . Complete study of K t4 is beyond the scope of this paper. In Ref. [19] we have presented 
a bosonized axial-vector current (Eq. (72) of Ref. [19]) which has been used to study many r 
mesonic decays. Theoretical results agree well with data. Using the axial-vector current and 
the vertices KaK*tt, K*Ktt, KApk, pKK, and pnn presented in Ref. [19] in the chiral limit, 
the form factors of the matrix element of the axial-vector current of K 14 

< ^ + (p+)^~(P-)\^\K + (k) >= — -{F(p + +p_) fl + G(p + -p-) fl + R{k-p + -p-) tl } (14) 

m K 

can be calculated. The calculation is similar to r —» IG-nv^Kpv done in Ref. [19]. Using 
Eqs.(40,43) of the form factor G of Ref.[18], it is derived 

n ( 13 ) 
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The same L 8 is obtained by calculating related term of the form factor iq + defined in 
Ref. [ 17 ]. 

From Eqs.( 12 ) the Li is obtained 


r 1 r 7 9 3 2 c. o 

L\ = -{9c-c 2 H--(1-) }. 

8 iy 2 8vr 2V g ’ 1 


It is learned from Ref.[10] that g 2 ~ 0(N C ), / 2 ~ 0 (N C ), and m 2 ~ 0 ( 1 ) in large N c 
expansion. Therefore, Li j2 ,3 ~ O(Nc)- The predictions of Li _ 3 are made at the tree level, 
therefore, the Eq.(12) actually means that 2 (Li + L 2 ) + L 3 ~ 0 ( 1 ) in large iVc expansion. 

According to Ref. 1 (b), the coefficients L 4 _ 8 are determined from the quark mass expan¬ 
sions of m 2 , m 2 K , m 2 , fx, and f v . Therefore, the quark mass term of the Lagrangian( 5 ) 
is needed to be taken into account in deriving the effective Lagrangian of mesons. The ex¬ 
pressions of the masses and decay constants of the pseudoscalars are found from the real 
part of the effective Lagrangian. LIsing the method presented in Ref. [ 10 ], in Euclidean space 
the real part of the effective Lagrangian of mesons with quark masses is written as 



H D n 00 1 1 

d D x - V -_-_ 

( 27 r) D n {p 2 + m 2 ) n 


Tr {(7 ■ d — ■ v + i'y ■ 075) (7 ■ d — i'y ■ v — i'y ■ aj 5 ) + 2 ip ■ (d — iv — ia ) 


+my • Du — i[y • v, M] + 7(7 • a, M }75 — m{uM + Mu) — M 2 } n . ( 17 ) 


where Dpi = dpi - i[v^u\ + f{a M ,u}75 and u = exp {-*75 [7*71-, + \ 0 K a + X 8 g 8 + 770]}. 
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Comparing with Eq.(ll) of Ref.[10], there are new terms in which the quark mass matrix M 
is involved. As done in Ref. [10], the mesons fields in Eq.(17) are needed to be normalized. 

The masses of the pseudoscalcar mesons are found from the coefficients of the terms 7777 , 
K a K a , r] 8 r]s in the effective Lgarngian(17). The masses of the pseudoscalars can be calculated 
to any order in quark masses. Up to the second order in quark masses we obtain 

4 i p 2 

ml± = -pj—{ —- <ip'ip> (m u + m d ) - —{m u + m d ) 2 }, 

F r± d 4 

4 1 p 2 

m l° = -jri-n <^> (m u + ™d) - ~p( m u + m d )}, 

j 71-0 O ^ 

4 i p 2 

m K+ = -FT -{-o <W> {m-u + m s )-— {m u + m s ) 2 }, 

Jk+ 4 4 

4 i p 2 

m K° = JT ~{ _ 3 < V’V’ > ( m d + rn s ) - — (m d + m s ) 2 }, 

4 1-1 F 2 

m 2 = > o ( m « + m d + 4m *) ~ + m d + (18) 

fri 3 3 6 

where F 2 is defined in Ref. [10] as 


t-i2 Nc 2 

t — ——m 

7T 4 


d A p 


('p 2 + m 2 ) 2 ’ 


< F’F’ > is the quark condensate of three flavors. In terms of a cut-off A the quark condensate 
is expressed as(using Eqs.(40) and (43) of Ref. [10]) 




7 ■ p — mu m 3 Nc r A 2 . , A 2 


p- — m z 


47 t 2 m 2 


+ 1 )}. 


m 


(19) 
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The universal coupling constant g is expressed as[10] 


g2 ~ 6^5 - 2^ { ' 0 » ( ^5 + 11 + AjVl “ 1} ' (20) 

m z 

The decay constants of the pseudoscalars f n , fx, arid f v are defined by normalizing the 
pseudoscalar fields as done in Refs.[10,11]- In this paper we calculate the contributions of 
the quark masses to these decay constants to 0(m q ). Using the effective Lagrangian(17), the 
decay constants can be calculated to any order in quark masses. As indicated in Refs.[10,11], 
there is mixing between the axial-vector field and corresponding pseudoscalar field. The 
mixing results in the shifting of the axial-vector field a p 

a /t -> —- —dfjP, (21) 

9a 9a 


where P is the corresponding pseudoscalar field, g a is the normalization constant of the o /t 
field, and c a is the mixing coefficient. Both g a and c a are determined in the chiral limit in 
Ref. [10]. For pion and cq fields it is obtained by eliminating the mixing between pion field 
and a i field 


1 F 2 F 2 3 

g 2 rn 2 2 8 m 2 in 2 


)2 m(m u + m d )}, 


( 22 ) 


where m a is the mass of the a\ meson which is determined from the Lagrangian(17) as 


9a m l = F 2 + g 2 m 2 p + 6 gl a m{m u + m d ), 


(23) 
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where g^ a is expressed as [10] 


& = 92(1 •* 2 ^)' (24) 

Up to the first order in quark masses, the decay constant f 2 is obtained from the La- 
grangian(17) 

fl = /. 2 o(l + f~ (25) 

where 

/S, = e 2 (i-^)[io] (26) 

and 

' = <* - 7)' 1 - 2 ^) - 1 + < *+ > ^ ~ 7 )(1 -;>• (27) 

In the same way, f] <+ , /J-o, and / 2 g are found 


fx+ ~ / 2 o(l + / 


m u + m s 


m 


}, 




m 


i2 _ ,2 r, , 

■bis — '631/1 1 + J J 


(28) 

(29) 

(30) 


Substituting Eqs.(25,28,29,30) to Eq.(18), to the second order in quark masses the masses 
of pion and kaon are obtained 

4 i ^2 f i 

ml± = -pr{-- <4>4>> ( m u + m d )-—(m u + m d ) 2 + 2- < i/ji/j > —(m u + m d )' 2 }, 

Jn 0 3 4 3 rn 
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4 , 1 


/ 


ml o = — (m u + m d ) - —{ml + mj) + - < ^ > —{m u + m d ) 2 }, 

£ O m 


f2 i- o 
J 7r0 ° 

4 , 1 


/ 


m K+ = -Fri~o <^> (m u + rn s ) -— (m u + rn s ) 2 + - < ifrij} > — (m u + m s ) 2 }, 

fio 3 4 3 m 

4 i p 2 f _ i 

m 2 K o = ~p 2 ~( q <^> {m d + m s ) - — (m d + m s ) 2 + - < ^ > — {m d + m s ) 2 }, 

< = W + m* + 4ro.) - ^(m 2 + m 2 + 4m 2 ) 

f - 11 

+ - < yy > ~ — {m u + m d + 4m s ) 2 } 

3 9 m 

41-1 f 1 F 2 1 

= <ip'ip> ~{m u + rn d + 4m s ) + (—- < ipip > ——) ~{m u + m d + 4m s ) 2 

JZo 3 3 m3 4 9 


2F 2 1 9 F 2 

——[~{m u + m d ) - m s ] 2 - — {m d - m u ) 2 }. 
y z iz 


(31) 


The dependence of the meson masses(31) on the quark masses are the same as the ones 
presented in Ref. (1(b)). The mass difference of charged pion and neutral pion is found from 
Eqs.(31) 

m 2 ± - m 2 o = (1 - —) _1 (m d - m u ) 2 . (32) 

9 

Comparing with Eqs.(31) with the Eqs. (10.7,10.8) of Ref.(1(b)), following coefficients of the 
ChPT are predicted 


3T7 + Eg 


U — 0 , 

L 5 

L 8 = 
F 2 

-MB 2 ’ 


U — 0 , 

(33) 

flvf 

8 mB 0 ’ 

(34) 

F 2 

(35) 

MB 2 ’ 

L 7 = 0, 

(36) 


13 



where 


4 1- 

B o = 7r( _ o) < W> > ■ 


/, 


7r0 


Using Eq.(37), L 5 and L 8 are written as 


where 


U = 


1 . 2c. 2c, 9/ 1 

(1 --){(! --) 2 ( 1 - 


32 Q 


9 


9 


27rV 


7-(i--) + 4q(i-7)}. 


9 


7T 


L 8 — 


1536^ 2 g 2 




1 nl 7 , 

Q = -3^4 <<P<P> 


(37) 


(38) 

(39) 


(40) 


The Eqs.(19,20) show that Q is a function of the universal coupling constant g only. The 
numerical value of Q is 4.54. 

Both L 5 and L 8 are at O(Nc)- I/ 4 , L 6 , and L 7 are from loop diagrams of mesons and are 
at 0(1) in large Nc expansion. 

By inputting the values of m 2 + , m 2 K+ , m 2 K0 , and / 2 (f n = 0.184GeU) into the Eqs. (18,25) 
we are able to determine the values of the three current quark masses and F 2 . The universal 
coupling constant g = 0.39 is chosen to fit the decay rate of p —> e _ e + . The results are 

F 2 = 0.0537GeU 2 , (41) 

m u = 0.91 MeV, rn d = 2.15 MeV, m s = 52.2 MeV. (42) 


The values of the quark masses are less than their traditional values[20]. As pointed out in 
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Ref. [20], the absolute values of rn q are model dependent. The values of m q determined by- 
lattice gauge calculation are smaller too[21]. Using Eqs.(41,42), the f K is calculated to be 

f K = 0.215GeV, 

Ik = 1.17/,. 

The experimental data is fx = 1.22/,. The deviation of the theoretical result from the 
experiment is about 4%. To determine m v and f v the mixing between r}$ and rjo is needed 
to be taken into account. From the effective Lagrangian(17) the terms related to the masses 
of rj and r/' are derived 

1 1-1 F 2 

= — 23 < 00 > 3( m « + m d + 4 '"U) ~ + m 2 d + 4 ™s)} 

11-2 F 2 

2 ( 3 < 00 > 3 K + - — (m 2 + m^ + m 2 ) + m’} 

^2 - V^F 2 

-^o^s{— < -0-0 > (m u + m d - 2m s ) -—— (m 2 + m d - 2m 2 )}, (43) 

y JLz 

where m s is a parameter resulted by the U( 1) problem [ 22 ], in Eq.(43) r/o and r/g are not 
normalized yet. The mixing is introduced as 

rj 8 = r/cos# — rj'sin9 , r/o = rjsind + rj'cosd. (44) 

Eliminating the crossing term of r/r/ and inputting m r; ' to determine m g , we obtain 

4 1 - F 2 

m l = -^{cos 2 ^-- < 00 > (m u + m d + 4m s ) - — (m 2 + m^ + 4m 2 )]} 

Jn 9 6 
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2 F 2 . 0 0 

+sin 2 9[-- < i/ji/j > (mn u + m d + m s ) -— (m 2 + m 2 d + m 2 ) + rriz] 

y o 

a/2 1 - F 2 

+^-sin29[- < -0-0 > (m u + m d - 2 m s ) + — (m 2 + m 2 d - 2m 2 )]}, 

f 2 = F 2 ( 1 - —){1 + [(1 - — )(1 - 1 ) - l](sin 2 a—(m u + m d ) + cos 2 a—m s )} 

1 g g 2 n z g z m m 

1 2c c 4 i 

H—y(l- X 1 -)^72 _ { cos 26 , ( _ q) < W > (m u + m d + 4m s ) 

7 T g g 3// 0 3 

+sin 2 9 [— < > 2(m u + + m s ) + 3 mz\ 

3 y 

y/2 

H— —sin29 < 00 > (m u + m d — 2m s )}, (45) 

o 

where a = 9 + 35.26°(sm35.26° = 4^). The numerical results are 

9 = 7.84°, = 0.569GeV, m® xp = 0.548GeV, (46) 

fn = 0.24 Gey = 1.3/,. (47) 

The theoretical value of m v deviates from the experiment by 3.8%. The values of the mixing 
angle and f v can be tested by the decay g —> 2y. Because the quark mass correction has 
been taken into account in calculating f v , the effect of quark masses has to be included in 
calculating the decay amplitude. This study is beyond the scope of this paper. 

According to Refs.[1(b),3], L 9 and Li 0 are determined by < r 2 >, and the amplitudes of 
pion radiative decay, n~ —> e~^v 


L - II <r>> 


La — 


R 


32vr 2 F v ’ 


(48) 

(49) 
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( 50 ) 


1 F a 

~ 32 F 2 !^ ~ 


where R, F x , and F A are r a , hy ,and Ha of Ref. [3] respectively. In Ref. [19] we have derived 
the expression of pion radius as 


6 


< > =- \- 

^ 7 T 9 I o />9 


{ (l_^)2_ 47r 2 c 2 } 


' ' 9 

which agrees with the data very well. Using Eqs. (48,51), it is predicted[23] 


u = \ cg + 16^ {(1 " 7 )2 “ 4A2} ’ 


(51) 


(52) 


On the other hand, L 9 is determined by Eq.(49) too. The three form factors of the decay 
amplitude of n~ —> e _ 7 z/ are presented in our paper[19] 


F a = 


g 2 m, m 2 p 2c 
K - + 


1 m, m 2 
2v / 27r 2 / w m^ 1 

2c, 


Vz fn ml 1 g 1 tt 2 # 2 ^ g 2 ttV 

Using the mass formula of the cq meson(in the chiral limit) [10] 


pV _ 

2V2F 2 U : 

(53) 

>-2V 9 4-‘. 

(54) 

P 1 + V2cg^4[24]. 

J 7 V 

(55) 


1 , 2 F 2 


(1 " 2 vy )m; " + m <” 


it is obtained 


m P (1 1 * 

m 2 2 ir 2 g 2 ' 


-i 


= 1 - 


2c 


(56) 
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Substituting Eq.(56) into R(55), it is derived 


R = 


3V2 


m n f n < r 2 > 7 


(57) 


This is the expression obtained by applying PCAC to this process[25]. Therefore, the coef¬ 
ficient £9 derived from Eq.(49) is the same as the one obtained from Eq.(48). 

Using Eq.(56), the ratio |— is written as 

F A ^ 2c, 


— = (1 - — ) 2 . 
F v 1 g’ 


(58) 


The coefficient L w is found from Eq.(50) 


r 1 1 2 1 

£.0 - -ys + - 5^5 


O’ 


(59) 


Both £9 and £10 are at O(Nc). 

In the expressions of the coefficicnts(13,15,16,38,39,52,59) there are two parameters: g 
and c which have been determined from p —> e _ e + and the ratio Al. respectively. The 
coefficients of the ChPT are completely determined by this effective chiral theory. The 
numerical values of the coefficients of ChPT are presented in table 3. 

To summarize the results. The orders of the coefficients of the ChPT in large Nq expan¬ 
sion are predicted as £1,2,3,5,8,9,10 ~ 0(Nc) and £ 4,6 are 0(1). These predictions agree with 
Ref.[l(b)]. It is also predicted that £7 ~ 0(1) in this paper. The numerical value of £7 de¬ 
termined in Ref. [1(b)] (see Table 1) is smaller. It is interesting to notice that £i + |£3 = —1.4 
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Tabic 3: Predictions of the Values of the coefficients 


10% 

10% 

10% 

10% 

10% 

10% 

10% 

10% 

10% 

10 3 Lio 

3.0 

1.4 

l 

00 

bo 

0 

4.77 

0 

0 

-0.079 

8.3 

-7.1 


which is compatible with the value shown in Table 2. It is necessary to point out that the 
expressions of the coefficients L 4)5 j6 ) 7 j8 are derived from Eq.(31), however, Eq.(18) is used to 
fit the masses of the pseudoscalar mesons. It is found in this paper that this way leads to 
consistent determination of the quark masses. There are no new parameters in predicting 
the coefficients of the ChPT. The quark masses determined in this paper are closer to the 
lower bound obtained by lattice gauge calculation. Theoretical values of f K and agree 
with data within about 4%. 

The author wishes to thank B.Holstein for discussion. This research was partially sup¬ 
ported by DOE Grant No. DE-91ER75661. 
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